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We present evidence for Mott quantum criticality in an anisotropic two-dimensional system of 
coupled Hubbard chains at half-filling. In this scenario emerging from variational cluster approxi¬ 
mation and cluster dynamical mean-field theory, the interchain hopping t±_ acts as a control param¬ 
eter driving the second-order critical end point Tc of the metal-insulator transition down to zero at 
Pf/t ~ 0.2. Below tf, the volume of the hole and electron Fermi pockets of a compensated metal 
vanishes continuously at the Mott transition. Above t\, the volume reduction of the pockets is cut 
off by a first-order transition. We discuss the relevance of our hndings to a putative quantum critical 
point in layered organic conductors, whose location remains elusive so far. 

PACS numbers: 71.30.+h, 71.10.Pm, 71.10.Fd, 71.27.-|-a 


A subject of strong current interest in condensed mat¬ 
ter physics is the metal-insulator transition (MIT) [1] 
with a low critical end point Tc at which the Mott transi¬ 
tion ceases to be first order [2-4] . The nature of this crit¬ 
ical end point and its universality class is a long-standing 
issue. From general considerations, one expects it to be¬ 
long to the Ising universality class [5, 6], similar to the 
liquid-gas transition, with the double occupancy playing 
the role of a scalar order parameter of the transition. 
A canonical example is three-dimensional (3D) Cr-doped 
V2O3 where critical exponents extracted from electrical 
conductivity measurements very close to the critical end 
point Tc ~ 450K are consistent with the universality class 
of the 3D Ising model [7]. In contrast, similar experi¬ 
ments on layered K-type charge-transfer salts with sig- 
nihcantly lower ~ 40K have indicated unconventional 
Mott criticality [8]. They stimulated subsequent exper¬ 
imental studies either objecting the existence of uncon¬ 
ventional behavior [9] or supporting it [10]. Theoretical 
scenarios of the two-dimensional (2D) Mott transition are 
also controversial, ranging from ordinary Ising universal¬ 
ity [11-13] to unconventional critical exponents [14]. 

Recently, the question of the nature of the 2D MIT 
transition has been raised again as new experiments on 
K-type and palladium dithiolene organic conductors sup¬ 
port either unconventional criticality [2] or 2D Ising crit¬ 
icality [3], respectively. As the conductors studied in 
Ref. 2 possess low-T ground states with various broken 
symmetries, the unconventional Mott criticality seems 
to be generic and unrelated to the proximity to sym¬ 
metry broken states. Instead, the fact that the critical 
end point Tc is relatively low suggests quantum effects 
to become important and necessitates a physical picture 
contrasting the one building on classical phase transi¬ 
tions [15-17]. Furthermore, a possible support for the 


2D Mott quantum criticality comes from the dynamical 
mean-held theory (DMFT) [18, 19], which reveals unex¬ 
pected scaling behavior of the resistivity curves in the 
high-T crossover region T Tc- A stringent test of 
the link between this scaling behavior and the quantum 
criticality appears, however, to be impossible since the 
latter is masked in the half-hlled 2D Hubbard model by 
the low-T coexistence dome [20-23]. Moreover, various 
numerical studies hnd that Tc remains hnite also in the 
presence of lattice frustration [24-30]. Finally, while the 
effective suppression of Tc can be achieved with disorder, 
it requires the proper treatment of Anderson localization 
effects [31-33]. 

In this Letter, we propose a different route to account 
for a low critical end point T^ of the MIT. Considering 
the fact that quantum fluctuations are enhanced in low¬ 
dimensional systems with spatial anisotropy, we investi¬ 
gate, using two complementary state-of-the-art numerical 
techniques, the effect of anisotropic hopping amplitudes 
and try to locate the putative quantum critical point at 
T = 0 in the phase diagram. 

Model and methods. We study the frustrated Hub¬ 
bard model on a square lattice with an anisotropic hop¬ 
ping at half-filling: 

ij,a i i,a 

with chemical potential /r and Coulomb repulsion U. The 
hopping tij is t along the chains and t± between the 
chains. By tuning the ratio t±/t from 0 to 1, we bridge 
the limit of uncoupled one-dimensional (ID) Hubbard 
chains {t± — 0) and the isotropic 2D lattice {t±/t = 1). 
In order to remove the perfect nesting antiferromagnetic 
(AFM) instability of the Fermi surface (FS) [34] also in 
the 2D regime, which would lead to an insulating state 
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FIG. 1. Metal-insulator phase diagram of the half-filled Hub¬ 
bard model (1) as obtained by the zero-temperature VGA and 
CDMFT at T = t/40. (Top inset) The combined VGA and 
GDMFT estimate for the critical temperature Tc terminating 
the first-order MIT; Tc is driven down to zero at ~ 0.2, 
thus providing evidence for the quantum critical nature of the 
MIT. (Bottom insets) FS topology close to the critical inter¬ 
action Uc in different regions of the phase diagram indicated 
by arrows. 


at any finite value of U [35, 36], we add geometrical frus¬ 
tration via next-nearest-neighbor hopping t' = 

The results are obtained by two complementary 
quantum cluster techniques [37], which can both be 
described within the framework of self-energy func¬ 
tional theory [38]. In the cluster extension of DMFT 
(CDMFT) [39], a cluster of Nf. interacting impurities is 
dynamically coupled to an effective bath. The impu¬ 
rity problem is solved using the quantum Monte Carlo 
(QMC) Hirsch-Fye solver and coupling to the bath is de¬ 
termined self-consistently. To also make the study com¬ 
putationally manageable down to the lowest temperature 
T = f/40 in the 2D regime, where the sign problem ham¬ 
pers the usage of the QMC solver on larger clusters, we 
use a 2 X 2 plaquette cluster. The 2x2 cluster is a minimal 
unit cell which allows one to capture both the ID umk- 
lapp scattering process opening a gap in the half-filled 
band [40] and short-range 2D AFM spin fluctuations. To 
trace the Mott transition at zero temperature, we use the 
variational cluster approximation (VCA) [41, 42] with a 
2 x2 cluster and one additional bath site per correlated 
site as a reference system [22], i.e., an effective eight-site 
cluster. In VCA, the grand potential fl is approximated 
by the self-energy functional (SEF) at its saddle point. 
As variational parameters, we choose the hybridization 
V between correlated and bath sites and chemical poten¬ 
tials of the reference system ^l' and the lattice system /r, 
respectively [43]. 

Phase diagram. Our main results are summarized in 
the ground-state phase diagram in Fig. 1. It shows our es¬ 
timate of the critical interaction strength Uc at which the 



FIG. 2. VGA self-energy functional fl in the proximity of 
MIT as a function of hybridization V/t for (a) t±/t — 0.2 and 
(b) t±/t = 0.3. In (b) stable minima are indicated by thick 
arrows; thinner ones mark unstable solutions. 


system undergoes a transition between Mott insulating 
and metallic phases in the full range between the ID and 
2D regimes. In agreement with the exact Bethe ansatz 
solution [44] and bosonization approach [45], VCA yields 
the Mott phase for any [/ > 0 in the ID limit [46]. As 
shown in Fig. 1, this changes dramatically upon coupling 
the chains: single-particle hopping t± shifts the criti¬ 
cal interaction Uc towards a finite value, thus enabling 
the interaction-driven MIT. Initially, Uc increases steeply 
with t± and then continues to grow nearly linearly, as ex¬ 
pected for the MIT controlled by the ratio of Coulomb 
interaction to kinetic energy gain. For t±/t > 0.2, the 
MIT line is found to be first-order consistent with for¬ 
mer studies of the frustrated 2D Hubbard model [24- 
30]. In contrast, in the strongly anisotropic case with 
tj_/t < 0.2, it marks a smooth metal-insulator crossover 
down to T = 0. This is supported by a systematic reduc¬ 
tion of the critical end point Tc identified within CDMFT 
(see the inset of Fig. 1). All of these aspects consistently 
suggest that t_L is a control parameter which tunes the 
nature of the Mott transition from strong first order in 
the 2D limit to continuous at ~ 0.2. We complement 
the phase diagram by showing the change of the FS topol¬ 
ogy when tuning t_L for values of the interaction close to 
Uc- Two main features come into play: (i) finite t_L leads 
to a warping of the ID FS and, in the presence of interac¬ 
tions, to the formation of hole and electron Fermi pockets 
of a higher-dimensional compensated metal [47], and (ii) 
for values t_L/t > 0.7 the compensated metal structure of 
the FS disappears going over to a conventional large FS 
which coincides with the topological Lifshitz transition 
of the noninteracting FS. 

Obtaining the phase diagram. We now describe the 
numerical results which lead us to the above phase 
diagram. VCA provides the possibility of identifying 
and tracing competing phases by analyzing the self¬ 
energy functional fI(/4, p,' ,V). For the interchain coupling 
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FIG. 3. (a) VGA ground-state energy Eq as a function of 
Goulomb repulsion U/t. Filled (dashed) lines indicate metal¬ 
lic (insulating) solutions for t^/t = 0.2,0.3,0.5,0.7, and 0.9 
(from left to right), (b) VGA double occupancy d across 
the MIT at T = 0; symbols stand for GDMFT results 
at T = t/40. (c) Cluster spin susceptibility Xail) within 

GDMFT at T = t/40; the upper curves correspond to spin 
fluctuations at the AFM wave vector q = (tt, tt) and the lower 
ones to remnant ID fluctuations at q = (7r,0). 


tj_/t = 0.2, we cannot resolve two disjoined SEF minima 
and the value of V is thus expected to change continu¬ 
ously across the critical interaction Uc, see Fig. 2(a). In 
contrast, for t±/t > 0.3, the SFF has four saddle points, 
of which two correspond to stable phases close to the 
phase transition: one corresponding to the metallic, the 
other to the insulating solution. These stationary points 
of ,V) are maxima with respect to and /i', but 

minima with respect to hybridization strength V. The 
existence of two distinct minima in the SFF landscape 
shown in Fig. 2(b) results in a jump of hybridization V 
when tuning across Uc, and thus signals the Hrst-order 
nature of the MIT. 

Next, we focus on the ground-state energy Fq and the 
double occupancy d. The latter is obtained as the deriva¬ 
tive of the grand potential 12 with respect to Coulomb 
repulsion d = In the quasi-2D region we identify a 
clear kink in the ground-state energy Eq] see Fig. 3(a). It 
arises from a level crossing of the insulating and metallic 
solutions and gives rise to a jump in the double occu¬ 
pancy d at Uci as shown in Fig. 3(b). The latter exhibits 
hysteresis in the region with two solutions, as expected 
for the first-order transition. Although a weak kink in Eq 
is also resolved for intermediate values of tj_, both the co¬ 
existence region and the jump in the double occupancy 
shrink and vanish at ~ 0.2 [48]. The absence of a 
jump in d and a single minimum in the SFF yield strong 
evidence for the continuous nature of the MIT. A similar 
scenario emerges within a finite-temperature CDMFT: 
while a clear jump in d = ^ found in the 



FIG. 4. Double occupancy d as a function of interaction U/t 
obtained in GDMFT for (a) t±/t — 0.3 and (b) t±/t = 0.9. 
The low-T jump in d signaling the first-order MIT turns into 
a crossover above the critical end point E- 

quasi-2D regime, it gradually decreases when reducing t_L 
and finally converts into a crossover at = 0.2. It re¬ 
mains smooth down to our lowest temperature T = t/40; 
see Fig. 3(b). As shown in Fig. 3(c), the level crossing 
in the ground state is also reflected in the spin sector 
and produces a jump in the cluster spin susceptibility 
Xs{q) = wjo at the AFM 

wave vector q = (7r,7r). In contrast, no distinction be¬ 
tween the response in Xs{q) at g = (7r,0) and q = 
wave vectors at t±/t = 0.2 indicates that remnant ID ef¬ 
fects play an important role in the weakly coupled regime. 

We turn now to finite-temperature consequences of the 
continuous MIT seen at T = 0. The estimate of Tc at 
a given tj_ was obtained by monitoring d as a function 
of U/t at fixed T; see Fig. 4. The low-T jump in d sig¬ 
naling the first-order MIT remains up to T^ and turns 
into a smooth crossover above Tc- As shown in Fig. 4(a), 
for small t±jt = 0.3 a smooth behavior in d is already 
recovered at T = t/30. In contrast, for t±/t = 0.9, the 
jump converts into a crossover at much higher tempera¬ 
ture T = t/12. By repeating the above analysis for inter¬ 
mediate values of t± , we extracted the t± dependence of 
the critical temperature Tc (see the inset in Fig. 1) con¬ 
sistent with a continuous reduction of the jump in the 
double occupancy [48]. 

Spectral function. To elucidate the microscopic ori¬ 
gin of the continuous Mott transition for tj_/t < 0.2, we 
calculate the single-particle spectral function A{k,uj) = 
— ilmG'(A:, w), where G{k,oj) is the lattice Green’s func¬ 
tion. Figure 5 shows the evolution of A(k,u}) upon in¬ 
creasing the interaction U at fixed t±/t = 0.2. In agree¬ 
ment with random-phase approximation studies [50, 51], 
we find that the destruction of the FS starts at mo¬ 
menta k = (7r/2, ±7r/2), where the interchain hopping 
matrix elements vanish. As a result, a compensated 
metal structure of the FS is formed with elliptic electron 
and hole pockets around the k = (7r/2,0) and (7r/2,7r) 
points. A striking feature of the pockets is their sym¬ 
metric form contrasted with pockets found in coupled 
spinless fermionic chains [52]. We ascribe this symmetry 
to quasiparticle scattering off short-range ID spin fluc¬ 
tuations with q = (tTjO). On one hand, at intermediate 
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(a) Ult = 2.7 (b) C//t = 3.2 (c) U/t = 3.22 



FIG. 5. Evolution of the FS with electron and hole pockets 
(see the text) for t±/t = 0.2 when approaching Uc/t ~ 3.22 
from below: (a) U/t = 2.7, (b) U/t = 3.2, and (c) U/t = 3.22. 
Red solid lines show the noninteracting dispersion. 


interaction strengths, the main part of the FS carrying 
most of the quasiparticle weight follows closely the non¬ 
interacting FS. On the other hand, the pockets shrink in 
size and become very shallow close to Uc, see Figs. 6(a)- 
6 (c). The continuous vanishing of the volume of Fermi 
pockets at C4 implies the second-order nature of the MIT. 
Since the inverse width of the hole or electron pockets de¬ 
fines a characteristic length scale, one should be able to 
extract the correlation length exponent, v, from a careful 
study of the critical behavior of the volume of the pocket 
as one approaches [1, 48]. In contrast, the volume re¬ 
duction of the pockets is cut off by a first-order transition 
for larger cf. Figs. 6(d)- 6(f) with t±/t = 0.5. 

Discussion and outlook. Let us relate our find¬ 
ings to recent experiments on the organic conduc¬ 
tors with a half-filled band [2]. Both k-(BEDT- 
TTF) 2 Cu 2 (CN )3 and EtMe3Sb[Pd(dmit)2]2 (BEDT- 
TTE=bis(ethylenedithio)-tetrathiafulvalene, dmit=l,3- 
dithiole-2-thione-4,5-dithiolate, Me=CH 3 , Et=C 2 H 5 ) are 
thought to be layered systems with Hiickel parameters 
close to an equilateral triangular lattice [53]. Instead, 
careful ab initio calculations for the latter show an ap¬ 
preciable ID anisotropy with the ratio of interchain to 
intrachain transfer around 0.82 [54]. 

We took this asymmetry into consideration: using 
VGA at T = 0 and CDMET at finite T we have found 
strong evidence for Mott quantum criticality in coupled 
Hubbard chains at half-filling. In this scenario, the in¬ 
terchain hopping t_L acts as control parameter driving 
the second-order critical end point of the interaction- 
driven MIT down to zero in the presence of strong 
anisotropy. At a threshold value of f^/t ~ 0.2, the vol¬ 
ume of Fermi pockets shrinks to zero. The resulting MIT 
is continuous without a detectable jump in the double 
occupancy or a visible coexistence region in the SEF. In 
contrast, the volume reduction of the pockets is only par¬ 
tial at larger t±: the jump in the double occupancy and 
the existence of two distinct degenerate minima in the 


SEF are consistent with a first-order transition. 

The continuous MIT at T = 0 offers a possibility for 
understanding the scaling behavior of resistivity curves 
in the high-T crossover region T ^ Tc usually attributed 
to (hidden) 2D Mott quantum criticality [18, 19]. It is 
an interesting question as to whether the quantum crit¬ 
ical behavior emerges also in coupled spinless fermionic 
chains displaying similar FS breakup into Fermi pock¬ 
ets [52]. 

While the 2x2 plaquette cluster used is known to 
overestimate the singlet formation [37], we expect the 
unveiled quantum critical behavior to be robust. Indeed, 
former CDMFT studies on larger clusters of up to 16 sites 
have provided evidence for a continuous dimensional- 
crossover-driven MIT down to the lowest accessible tem¬ 
peratures [55]. We believe that this scenario is not re¬ 
stricted to quantum cluster descriptions of the system 
but should also emerge in lattice simulations, provided 
that the range of AFM spin fluctuations is reduced, e.g., 
by geometrical frustration or disorder [56, 57]. This leads, 
however, to a severe sign problem which renders lattice 
QMC simulations very expensive [58]. In this respect, a 
promising route avoiding the main shortcomings of QMC 
is offered by tensor network methods [59] adapted re¬ 
cently to fermionic systems [60, 61]. Our results provide 
a novel axis in the phase diagram along which Tc can 
be tuned to zero. It remains to be verified if this quan- 
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FIG. 6. Low-energy part of the single-particle spectral func¬ 
tion A{k,Lj -I- ir]) for t±/t = 0.2 (top panels) and t±/t = 0.5 
(bottom panels) obtained within VGA at T = Q, rj = 0.05. 
The FS pockets are found (a),(d) for U < Uc and (b),(e) at 
U < Uc, in the insulator at (c),(f) U > Uc, their disappear¬ 
ance signifies the vanishing of the FS, and hence the MIT. 
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turn critical behavior can explain fingerprints of the un¬ 
conventional Mott criticality observed in layered organic 
conductors. 
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Supplemental Material for: Mott Quantum Criticality 
in the Anisotropic 2D Hubbard Model 

Phase diagram Fermi surface pockets 


In VGA we discriminate metal against insulator by 
looking at the electron filling (n) as a function of chemi¬ 
cal potential /r. It shows a plateau at half filling for large 
Coulomb repulsion U, which corresponds to a vanishing 
charge compressibility k = — and indicates an insu¬ 
lator, whereas for a metal no pfateau occurs. 


Double occupancy 


Within CDMFT with the QMC solver the model sys¬ 
tem can be analyzed at hnite-temperatures T. This al¬ 
lows one to estimate the t_L dependence of the critical 
endpoint temperature T^. In Fig. IS we show the double 
occupancy d as a function of Coulomb repulsion U/t for 
intermediate values of the interchain hopping t^/t = 0.5 
and = 0.7 thus complementing Fig. 4 in our Letter. 

As there is no jump in d within CDMFT for small inter¬ 
chain couplings t±/t < 0.2 down to the lowest accessible 
temperature T = t/40, VGA is used to investigate the 
system at T = 0. Figure 2S shows the size of the jump 
in d as a function of interchain coupling. For interme¬ 
diate and large t±, the metallic and insulating solutions 
both exist in close vicinity of the transition and the jump 
can be read off immediately. In the case of small t±, it 
is necessary to fit double occupancies of the metallic and 
insulating solutions to read off the jump at the transition 
(see inset in Fig. 2S). Although an exponentially small 
jump for t±/t < 0.2 can still not be ruled out in this 
way, a fit of the jump sizes for t±/t > 0.2 agrees with a 
vanishing jump size for small t±. 



Figure 3S shows the FS for different values of tj_/t close 
to the critical value Uc- When comparing to the tight- 
binding limit, the main effect of the interaction is the 
formation of additional shadow-like bands, and the cre¬ 
ation of gaps in the originally continuous structure. Both 
effects lead to the formation of the compensated metal 
structure of the FS with hole and electron pockets in the 
region tj_/t < 0.7. In the case of a continuous transition 
ioi t±/t < 0.2, the volume of the hole and electron Fermi 
pockets shrinks to zero at the critical interaction strength 
Uc as found in both VGA and CDMFT, see Figs. 4S(a,d). 
For larger t± the pockets retain a finite volume up to 
Uc, see Figs. 4S(b,e) with t±/t = 0.3 and 4S(c,f) with 
tj_/t = 0.5. In the level crossing scenario, single-particle 
spectral weight is substantially redistributed when cross¬ 
ing from metallic to insulating solutions. This is seen as 
a jump in the single-particle gap as well as other observ¬ 
ables indicating a first-order transition. 

To quantify the interaction-driven renormalization of 
the FS warping for t±/t = 0.2, we plot in Fig. 5S the 
width of the electron pocket as a function of Ujt. For 
U ^ Uc, the main part of the FS carrying most of the 
quasiparticle weight closely follows the tight-binding dis¬ 
persion. Increasing the interaction strength results in a 
measurable deviation from the tight-binding dispersion 
and the continuous metal to insulator quantum phase 
transition corresponds to the shrinkage of the volume of 
Fermi pockets on approaching the critical point Uc- Since 
the inverse width of the hole or electron pockets defines 
a characteristic length scale, one should be able to ex- 
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FIG. IS. (Color online) Double occupancy d as a function 
of Coulomb repulsion Ujt obtained within CDMFT for: (a) 
= 0.5 and (b) tj_/t = 0.7. The low-T jump in d signaling 
the first-order MIT turns into a crossover above the critical 
endpoint Tc. 


FIG. 2S. (Color online) Size of the jump in double occupancy 
in VGA at the metal-insulator transition at T = 0. The 
dashed line is a rough fit and only meant as a guide to the 
eye. For t±/t < 0.2 we cannot identify a jump in d anymore 
and the transition seems to be continuous. 
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FIG. 3S. Fermi surface topology for several values of the inter¬ 
chain coupling t±/t close to the critical interaction strength 
Uc- For tj_/t > 0.7 the compensated metal structure of the 
FS with hole and electron pockets disappears going over to a 
conventional large Fermi surface. 
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tract the correlation length exponent v [2]. Thus, the 
analysis of the critical behavior of the volume of Fermi 
pockets provides an opportunity to address the univer¬ 
sality class of the MIT along the quantum critical line 
below = 0.2 which is left for future. 

The change in the nature (second- vs first-order) of 
the MIT with increasing interchain coupling is also seen 
in the t_L dependence of the pocket width as a kink at 
tj_/t = 0.2, see Fig. 6S. It marks the crossover between 
the regime where the volume of hole and electron Fermi 


FIG. 5S. (Golor online) Size of the electron pocket for t±/t = 
0.2 as a function of U/t. Length of the pocket Ly has been 
divided by a factor of two. The tight-binding limit [/ —>■ 0 is 
indicated by the dashed lines. 


pockets vanishes continuously at the second-order Mott 
transition and the region where the volume reduction of 
the pockets is cut off by a first-order transition. 
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FIG. 4S. Low-frequency part of the single-particle spectral 
function A{k,LjJ + iri) in the metallic phase close to the critical 
interaction Uc for different t± obtained within VGA at T = 0 
(top) and GDMFT at T = t/40 (bottom): (a,d) t±/t = 0.2; 
(b,e) t±/t = 0.3, and (c,f) t±/t = 0.5. In VGA small broaden¬ 
ing rj = 0.05 was used; stochastic analytic continuation [1] of 
the imaginary-time QMG data was applied within GDMFT. 
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FIG. 6S. (Golor online) Width of the electron pocket as a 
function of t±/t close to Uc- The location of the kink agrees 
with the change in the nature of the MIT. 









































































